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IIYDROBLASTIC EBBECTS IN THE AORTA BIFURCATION ZONE 


A. S. VolVmlr, M. S. Qerahtoyn and B. A, Purinya 

Ain Pox-*ce Enc;lnoerlnR Aoademyj mobcow) 

Latvian SSR Academy of ScienceB 
Institubo of MoehanicB of aolymora* Riqa 

Accoi-’dinn to data from aatopalGs ond cardiovascular surEer>y, /Id^l 
the procoQS of sclerotic plaque formation and embrittlement in persona 
of all ar.cM , but chiefly in the elderly, frequently occurs where the 
major vessels ramify, especially at the aortic bifurcation. Serious 
diseases resultiuK in inadequate blood supply to the Igrs are asso- 
ciated with those processes. 


Analysis of the hydroelastic aspects of such phenomena would 
lead to a better under standinR of their orif.lno. One obvious cause 
of atherosclerotic plaques, most often encountered at the edges of 
the bifurcation, is the foimiation of closed eddies where other major 
vessels branch off from the aorta. Combined solution of hydrodynamic 
and elasticity equations is required for analysis of the mechanical 
behavior of blood and vessels at the bifurcation. Our work [1] has 
shown that, in general, a blood vessel should be treated as a multi- 
layered envelope, one layer of which — the muscular — through active 
contraction exerts a substantial influence upon the ’nitial (for each 
given hemodynamic problem) condition of -the vessel-envelope, 

Structurally, the aorta is an elastic-type vessel. The inner 
coat of the aorta (the intima) is distinguished by its complex 
structure and, in comparison with the intima of smaller-caliber 
vessels, its relative thickness. The middle coat of the aorta is 
made up of docens of alternating layers of elastic and smooth muscle 
tissue, and its external coat is comparatively poorly developed. 
Bearing in mind that the intima is considerably less rigid than the 
other coats, while the middle coat, because of its large number of 
layers, is essentially homogeneous, we can use a single-layer envelope 
as a model of the aorta. Similarly [P], we shall write the equation 

^♦Numbers in the margin indicate pagination in the foreign text. 


1 



of motion for it as; 
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where ,v, y, anc'. ^ are the coordinates for the median surface of the /l6g 

envelope. Displacements in the x, y, and a directions are correspondingly 

designated as u, v, and ?*), median surface stress as N^. and T, and 

^ y 

the bending and torque moments as M , M , and H. The density and thick- 

X y 

ness of the envelope material are p and h’, k and k are the median 

X y 

surface curvatures. 


Now let's deal with the blood flow equations. A number of models 
are known to describe blood rheology in the cardiovascular system [33. 
When analyzing blood flow in the major vessels, we may consider the 
blood a viscous Newtonian fluid whose movement is described by the 
Navier-Stokes equations 
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(5) 


and the continuity equation 

dp/ 0(piVr) ■ I d(p/oo) d(p/t>»;l . d(pyM 
d/ dr r ■ dO *“ *“ ox dr 

where a, r, and 6 are fluid particle coordinates 5 and Ug are 

its velocity projections j v « u/p^ is the kinematic-viscosity coeffic- 
ient j is the fluid density. The equatilon of state p "w p(p^) should 
be added to these equations. 


Boundary conditions for the fluid at the deformed walls of the 
vessel are: 


Vr» 


doy 

'"It 




du 

Ot 


ao* 


jdt 

dF 


( 6 ) 


The forces exerted by the flowing, blood on the walls of the aorta 
are determined using the relationships 
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As can be seen in equations (1) - (3), we shall disregard the force p . 

y 

The stationary blood flow through the bifurcation during diastole 
may be taken as the original state of the system. Combined Integration 
of the system of equations for the envelope (1) - (3) and the hydro- 
dynamic equations (H) and (5) under the corresponding envelope and 
fluid resistance conditions (6) and (7) should be performed, giving 
the changes in pressure at the entrance to the area of the aorta under 
consideration at the bifurcation and the blood flow during systole. 

These "entry" functions may be provided as based upon measurement data 
which are available in the literature. 

The difficulties associated with integrating these equations 
for a non-symmetrical area and in the event of a non-stationary process 
are obvious. They can only be surmounted by using numerical methods /166 
and a large digital computer. 
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As a first approximation, we can 
ccmsicler the problem of fluid flow in a 
similarly shaped smooth channel, which 
will provide a tentative description of 
blood flow throuRh the bifurcation. This 
problem is considered C53 with the assumption 
that the channel walls are rigid, Navier- 
-Stokes equations for an incompressible 
fluid are transformed into the eddy 
transfer equation C6]: 
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where is the flow function assoblated with the fluid velocity ratios 
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When velocity changes are uneven, a closed eddy gradually forms at the 
external angle of the bifurcation near the entrance to the large 
channel, but it substantially affects the primary flow sone as well 
(fig. 1). The most intense fluid particle rotation is found at the 
outside of the streamline, going away from the separation point. 


As the current develops in time, the streamline nearest the wall 
gradually diverges from it , while the separation point (designated 
in fig. 1 by 0) shifts upward somewhat along the flow. A low pressure 
area forms at the center of the closed eddy. Of extreme importance 
in this case are the high shift stresses in the area where the closed 
eddy forms. 

The sclerotic phenomena mentioned above are obviously associated 
with this condition. 


We should note that similar processes are observed in the ducts 
of large hydraulic systems, and they consist of salt deposition on the 
walls where the duct systems branch. 
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